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The Advanced Receiver (ARX II) currently being developed uses a Costas cross- 
over loop to acquire and track the phase of an incoming quadrature phase-shift- 
keyed (QPSK) signal . This article describes the performance of the QPSK lock 
detector to be implemented, taking into account the phase jitter in the tracking 
loop . Simulations are used to verify the results of the analysis. 


I. Introduction 

The Advanced Receiver II (ARX II) [1] is currently 
being developed to demodulate signals from deep space 
spacecraft. In addition to processing binary phase-shift- 
keyed (BPSK) signals, the ARX II will acquire and track 
quadrature phase-shift-keyed (QPSK) signals. The track- 
ing as well as the acquisition performance of several QPSK 
loops has been investigated [2], and it was determined that 
the Costas cross-over loop offered a “good” compromise 
between implementational complexity and relative perfor- 
mance. The other two candidates were the maximum a 
posteriori (MAP) estimation and the generalized Costas 
loops. This article describes a lock detector for the QPSK 
loops. The analysis is general in that it is applicable to 
all three QPSK loops, and the results are validated by 
computer simulations. The general QPSK carrier-tracking 
loop structure and the lock detector are both shown in 
Fig. 1. When this loop is implemented, an extra accu- 
mulator is present in front of the loop filter to reduce 
the loop update rate and improve the loop signal-to-noise 
ratio (SNR). The lock detector structure can be gener- 


alized to accommodate any multiple phase-shift-keying 
(MPSK) signal as described in [3]. 

II. Lock Detection Analysis 

The received QPSK signal can be modeled by 
r(t) = \fPoa{t) sin(u>o2 + 0) 

-F y/ P dK t) cos(ljo t + 0) + n(t) (1) 

where 

a (0 = ~ 


b(t) = ^ b k p(t - kT s ) 

fc 

are the independent in-phase and quadrature data streams 
with ajfc and bk the binary ±1 random data and p(t) the 
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non-return-to-zero (NRZ), or Manchester, pulse of dura- 
tion over T s seconds. The data power is given by Pq and 
the incoming phase by 0. The signal is received in the 
presence of noise that can be expressed as 


where z is the detector’s signal obtained from samples i /k 
Vk = ( r 2 ck - r] k ) 2 - (2 r ck r, k f (5b) 


n(<) = \Pi n c (t) cos(u>o t + 0) — V2n,(t) sin(wo< + 0) (2) 

where n c (t) and n 3 (t) are random processes with a two- 
sided noise spectral density of No/2. The input signal is 
mixed with the reference signals 


and 


( 4 \ 2VPd ( 4 , jj\ 

X c (t) - — — cos(u 0 t 4- 0) 
Nq 


X s {t) = sin (u> 0 t 4 * 0) 

No 


(3a) 


(3b) 


to produce, after integrating and dumping, the samples 


It is shown in Appendix C that the mean and the variance 
of the detector signal z can be expressed as 


= — 4MR 4 cos4(j> (6a) 


and 


2 M 


mg + 1977144 — 12 m <32 
1 M — 1 

4 - 4 A/ | — Cfcdk — Mm\ — 9 mi ~b 6/7/477?? 


k — 1 


(6b) 


r s k 


< 

=r ( 


(* + l)T. 


r(t)x 3 (t) dt 


1 1 • 1 Xo k * 1 

ak cos <p — bk sin 0 — cos 0 -= sin <p 

VR v R 


wdiere 0 denotes the carrier phase error; m, = E{r*} = 
E{r*} the ith moment; = E{rJ.r^} the ij th cross- 
moment; and cjt, djt are constants defined in Appendix C. 
(The overbar and E indicate expected values.) Note that 
at high-loop SNR (i.e., when a $ — ► 0), the mean and vari- 
ance reduce to 


r(k+i)T. 

r c k = / HOMO dt 

JkT, 


■ H ( 


Xjjt . X2k 

a jt sin 0 + bjc cos 0 — sin 0 4- cos0 

VR VR 


( 4 ) 


where R = ls ^ ie s y m bol SNR, Xu and xofc are 

two zero-mean, independent white Gaussian random se- 
quences with unity variance, and 0 = 0 — 0 is the phase 
estimation error. The reference signals of Eq. (3a) and 
Eq. (3b) need to be normalized by V Pd /N o for only the 
MAP estimation loop, which contains the “hyperbolic tan- 
gent” nonlinearity. For the other two loops, the amplitude 
can be a constant independent of V Pjy / No • The tracking 
performance of the three QPSK tracking loops has been 
derived elsewhere [3], and the results are summarized in 
Appendix A. 


The lock detector algorithm used in the ARX II is sum- 
marized in Appendix B for a general MPSK signal. In the 
specific case of a QPSK signal, the algorithm reduces to 

M 

z = X] yk < t ( 5a ) 

k = l 


fi z = - AMR 4 (7a) 

and 


64 A/ R 8 


'2 _9_ V2 _3_' 

R + R? + 7^ + VC 


(7b) 


from which the ideal detector’s SNR can be derived. At 
lower loop SNRs, the phase jitter might not be negligi- 
ble and might result in a degraded detector SNR. All the 
derivations have been summarized in Appendix C for ref- 
erence. 


III. Probability of Detection and of False 
Indication 

During carrier lock detection, each z sample is com- 
pared with a predefined threshold r, and the lock detector 
decides that the loop is in-lock when z exceeds r (i.e., 
z > t). It is possible that even when the loop is still not 
locked or when no signal is present, z will occasionally be 
larger than r. In this case, the lock detector will mistak- 
enly declare an in-lock condition. The probability of this 
event (probability of false indication) is 
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where p z o a nd a ?o are the mean an d variance of the lock 
detector signal in the out-of-lock state and erfc(x) is the 
complementary error function (erfc(x) = l-erf(x), where 
erf(x) is the error function). Equation (8) assumes that 2 
is a Gaussian random variable since z is the sum of many 
samples of equal variance. When the loop is not locked 
on the signal, the input can be modeled as pure noise to 
produce 


HzQ — 0 

<r 2 0 -192 MR 4 


( 9 ) 


On the other hand, when the loop is locked on the signal, 
the probability of detection is given by 


Pd = 



) 


dz 


— erfc 

2 


T- Vz \ 


( 10 ) 


where and <j 2 are given by Eq. (6). Alternatively, the 
threshold can be eliminated from Eq. (10) and the proba- 
bility of detection expressed as 


Pd = 2 erfc 


& z 0 


erfc“ 1 (2P / ) 


SNR Z 


( 11 ) 


where SNR Z - /iJ/o-J denotes the detector SNR. The 
phase jitter in the tracking loop degrades the detector s 
SNR by a factor D, 


D = 


SNR Z 

S N R z f ideal 


( 12 ) 


where SNR Z)id eai is the detector SNR if an infinite-loop 
SNR is assumed; i.e., there is no phase jitter (SNR Z (ideal) 
is computed by using the high-SNR approximations for fi z 
and a 2 given by Eq. (7)). In Appendix A, it is shown that 
at high-loop SNR the detector’s SNR reduces to 


S N Rz .ideal 


MR 

4(2-f| + ^+l^) 


(13) 


For a given M, loop SNR, p, (p = \/a\ where is given 
by Eq. (A-13)), and Pj, the detector’s SNR, has to be 
increased roughly by the factor 1/D in order to achieve 
the desired probability of detection P d . 


Given P d and Pj, we can solve for the number of re- 
quired detector samples M , namely, 



a y o erfc l (2Pj) 


T 2 

<7 y erfc _1 (2Pd) 


(14) 


The threshold r is obtained by solving Eq. (8) and setting 
it equal to 


r = y^2cr 2 0 erfc 1 (2 Pj) (15) 


where erfc -1 ( ) is the inverse complementary error func- 
tion. When the loop is in lock, it can be argued via the 
central-limit theorem that the random variable z is ap- 
proximately Gaussian, with mean and variance as obtained 
earlier. 


IV. Discussion and Numerical Results 

Computer simulations were used to validate the various 
assumptions made in the analysis. Figure 2 depicts the 
detection probability as a function of loop SNR when Pj = 
1CT 4 . It is clear that the detection probability degrades 
significantly at low-loop SNR and the degrad atio n itself is 
a function of the desired detection probability P d (which is 
computed assuming zero phase jitter). The degradation is 
better shown in Fig. 3, which depicts the loss in detector 
SNR as a function of the loop SNR. At 20-dB loop SNR, 
the degradation is about 1 dB, but it increases to 2.5 dB 
as the loop SNR is reduced to 15 dB. In both figures, the 
simulations and theory are in agreement for all three loops. 
Note that the degradation in detector SNR for BPSK was 
about 1 dB for 15-dB loop SNR [5]. 

When the detection probability was plotted as a func- 
tion of symbol SNR, different performances for the various 
loops were expected, as the loop SNRs are different when 
both the data rate and the loop bandwidth are fixed. Fig- 
ure 4 clearly depicts that effect, as the MAP estimation 
loop attains a higher detection probability than the Costas 
cross-over loop because of its higher loop SNR. Finally, 
Fig. 5 shows both the nominal and the actual (16-dB loop 
SNR) detection probabilities as a function of the number 
of samples M and the corresponding threshold. This kind 
of curve is very useful when one is designing the operat ing 
parameters of the lock detector. 
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Figure 6 depicts the probability of detection versus the 
detector’s SNR for several values of probability of false 
lock. This figure is useful in determining the required de- 
tector’s SNR for a desired detector’s performance. Finally, 
Fig. 7 plots the detector’s SNR versus E s /Nq for M — 300. 
The detector’s SNR for other values of M can be computed 
from the following relation: 

SNRz^m = SNRz'Zoo 4- 101og 10 d B (16) 

As a design example, suppose that a QPSK detector 
operates at Pj = 10 -4 and Pd = 0.99, and that the signal 
rate is 10,000 symbols per second with E s /Nq = 1.0 dB. 
Using Fig. 6, it can be seen that the detector’s SNR should 
be about 12 dB. But, Fig. 7 states that at E s /No = 1 dB, 
M = 300 and loop SNR = 16 dB (low-loop SNR case), 
the detector’s SNR will actually be about 5.0 dB. Equa- 


tion (16) is used to compute the number of samples, M, 
required to achieve 12 dB of the detector’s SNR, which 
results in M = 300 x 10( 12 ~ 5 O)/1 ° = 1504 = 0.15 sec, and 
Eq. (15) is used to compute the threshold r needed for 
Pj = 10~ 4 . Alternatively, Fig. (5) can be used to find the 
required threshold value, which turns out to be r = 1300, 
assuming that the output of the integrate-and-dump de- 
vices is scaled by 1.0/(E s /Nq). 


V. Conclusion 

This article analyzed the lock detector performance for 
the QPSK carrier-tracking loop in the ARX IT. The analy- 
sis is general and is applicable to other QPSIv loops. Both 
analysis and simulation were used to show that the effect 
of phase jitter on the detector SNR can be significant (as 
much as two decibels), especially at low-loop SNR. 
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PROBABILITY OF DETECTION {Pd) 



Fig. 1. Implementation of the QPSK loop with lock detector. 




Fig. 2. Probability of detection versus loop SNR. 


Fig. 3. Degradation In detector SNR versus loop SNR. 
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Fig. 4. Probability of detection versus £ s /Mq. (Theoretical per- 
formance assumes degradation due to phase jitter in the loop.) 
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NUMBER OF DETECTOR SAMPLES, M 

Fig. 5. Probability of detection versus number of detector samples. (The values of r are 
given with the assumption that the outputs of the integrate-and-dump devices are scaled by 
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Fig. 7. Detector SNR versus E s /N 0 . 




Appendix A 

Tracking Performance of Three QPSK Carrier Loops 


Depending on the loop, the samples r c and r, of Eq. (4) are processed to produce the error signal that forms the 
input to the loop filter. In general, that error signal can be expressed as 

= r ck f(r 3 k) - r, k f(r c k) (A' 1 ) 

where /(x) is some nonlinear odd function that defines the loop. Three possible functions will be considered. 

sgn(x) 


/(*) = { 


x° 

tanh(z) 


(A-2) 


In order to predict the performance of the loops, we need to compute the “squaring loss,” Si , which is defined by 


A ( E {^(°)}) 2 

E{6*(0)} 


(A-3) 


where e(0) and c r (0) are the error signal and its slope, respectively, evaluated at <j) = 0. Note that E{f 2 (0)} — N 2 equals 
the variance of the open-loop phase error process, From this, one obtains the equivalent single-sided noise spectral 
density, Nq , using the relation Nq = 2N 2 T S , where T s is the symbol time. The loop responds to the mean of the error 
signal, which is given by 


e(<j>) = E jf?(ajfc sin <j> + b k cos <j>)f (a k cos <f> — b k sin<£ + (-x lt cos <j> - x 2k sin <t>)/\fR} j j 


E i^R(a k cos <f> — b k sin <J>)f sin <j> + 6* cos<£ + (— x,t sin <j> + x 2k cos 


(A-4) 


where the expectation is taken first over the random data and then over the thermal noise. Using Eq. (A 4), the expected 
value of the slope of the error signal can be expressed as 


m = 




d<j> 


= RE^a k f [/?( ajb — X\k + b k f + X2le /Vk)]} 

+ R 2 E^b k r [*( <*k - Xik />/*)]( -bk - x 2 k /VR)) 
-R 2 E^a k f [f?(6* 4* X2k />/*)] (a k - Xik /VS)} 


(A-5) 


Averaging the above equation over the random data, one obtains 


c'(0) = 2 Rf[R - VRx] - 2 R 2 f'[R - VRx] 


(A-6) 
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where the overbar denotes expectation over the Gaussian random variable x with zero mean and unit variance. After 
some algebra, we find the expected value of e 2 (0), namely, 


E{e 2 (0)} = 2 R 


(R + 1)E {f 2 [R - y/Rx]} - (e {xf[R - y/Rx]} - y/~RE{f{R - VSe]}) ' 


(A-7) 


In order to proceed further, the non-linearity f(x) needs to be specified. For example, f(x) = sgn(x) (/'( x ) = 2<5(;r)) for 
the Costas cross-over loop, in which case we obtain 


E{sgn {R - y/Rx)} = erf (y/R/2) 

(A-8a) 

E{2 8(R - y/Rx)} = \f^ e ~ RP 

(A-8b) 

Efxsgn^- y/Rx)} = -\j\ e ~ R ' 2 

(A-8c) 

Esgn(R-y/Rx)} = 1 

(A-8d) 


which result in 




R + 1 - (y/R erf y/R / 2 + e-*/ 2 ) " 


(A-9) 


On the other hand, for the generalized Costas loop, f(x) = x 3 ( f f (x ) = 3r 2 ) and we have 


E {(R - y/Rx) 3 } = R 2 (R + 3) 

(A-lOa) 

E{3(f? - y/Rx) 2 } = 3 R(R + 1) 

(A-lOb) 

E{x(ft - y/Rx) 3 } = -R? I2 (R + 1) 

(A-lOc) 

E{/ 2 (x)} = R 3 (R 3 + 15 R 2 + 45 R + 15) 

(A-lOd) 

St — 1.9,6, 3 

^ OR ' ~r 7 ~r orJ 

(A- 11) 


In the MAP estimation loop, the nonlinearity is given by f(x ) = tanh(x) (/'(x) = sech 2 (x)), which prevents a closed-form 
expression, i.e., 




|tanh( J R — y/Rx) — Rsec\i 2 (R — y/Rx) 

2 

(R -f l)tanh 2 (/£ - y/Rx) - 

x tanh (R - y/Rx) — y/R t< 

anh (R — y/Rx) 


(A-12) 
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The variance of the phase error process in all three loops can be expressed as 


.2 = KBl _1_ = Bl 

* [e'(0)] 2 PcSl r ay RS L 


(A- 13) 


where p c = , Bl is the loop bandwidth, r, y is the symbol rate, and R = E,/Nq is the symbol SNR. Note that 

for R « 1 dB, Sl •‘C 1 . We should comment that the above variance u ^ is the variance relative to wo, which is used to 
compute the telemetry performance. The variance that determines acquisition and the cycle slip process is relative to 
4wo, and hence will be 16 times larger. 
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Appendix B 

Error and Detector Signals for MPSK 


Let r c be the “in-phase” and r, the “quadrature” outputs of the integrate-and-dump devices of a MPSK loop: Ref. [4] 
proposes the following algorithms for deriving the error signal and the detector’s signals for the various MPSK tracking 
loops. 


Derivation algorithms 


M 

Error signal 

Detector signal 

2 

e 2 = r c x r. 

d 2 = rl- r'i 

4 

e 4 = d 2 x C 2 

d-A — d 2 — 

8 

eg = d 4 x €4 

d s = d.\- e\ 

N 

CJV — d I± X € iv 

d\[ — d~ N — e^ 

2 2 
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Appendix C 

QPSK Lock-Detector Statistics 


The signal for the lock detector is derived using the following algorithm 

M 

z ~Yl yk> < T 

k 


where 


Vk = ( r ck - r,k) 2 - (2r c kr,k) 2 (C’ 1 ) 

In order to compute the mean and the variance of z y we need the second through the eighth moments of r c k and i s k- 
The odd moments are all zero (since the noise is Gaussian), and hence only the even moments need to be computed. 
To simplify the notation, the subscript k will be dropped from the next several equations. Taking the expected value of 
Eq. (4), the even moments are given by 



rri4 


= — rj = R 4 — icos4<^ 


m 6 



15 


15 


= ( o ( 5 - 3 c °s4 <f>) + —(3 - cos 4 <f>) + 


R 2 



(C-2) 

(C-3) 

(04) 


m$ - 


A .8 - r 8 = ft * 


16 


105 — 5cos4<£ — 84cos4<£ 


) 


/ 7 \ / 105 \ 420 105 

+ ( 2^(5 + 3cos 4^)+ 15(1 -cos 4^) jf — (3 + cos4^) + 3(1 - cos 4<^) j+ — + — 


(C-5) 


where 


r = 3 



2_ J_ 

R + R? 


At high-loop SNR, <f> — + 0 and m 4 , mg, mg reduce to 

RA ( ! + I + 7?) 

15 45 15 


R {l + ~R + W + R 3 


R a [l + — + 


28 210 420 


+ 


R R 2 R 3 


105 \ 

R 4 ) 


(C-6) 


(C-7) 

(C-8) 

(C-9) 
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respectively. Cross-moments for i ~ j are also needed and are given by 


m 22 = E{r 2 r 2 } = R 4 Q(1 + cos4<£) + J; + T 


18 


ni 44 = E {r 4 r 4 } — R s ( r 2 + -(1 -f cos8<^) — rcos4<?!> — — (1 — cos4</>) 


8 


(C-lOa) 

(C-lOb) 


m 62 = E{r, 6 r c 2 } = E{r®r 2 } = R 8 ( -5+ 6cos4^ + -^(25 - 9cos4<£) + ^(9 - cos 4^) + ^ + 


15 


60 15 


R ' 


2R' 


R 3 R 4 


(C-lOc) 


Note that as a# — > 0, m 22 — ► 77144 — ► 7 tz|, and rriQ 2 — ► 7riQm2. Now, the mean value and the variance of the detector’s 

signal can be computed. The mean is given by 


M 

Hz = E = M Py 

Jfc=l 


(C-ll) 


where 


Hy = r 4 + r 4 - 6r 2 r 2 


(C-12) 


Using the moments just obtained gives 


fly = ~4/? 4 cos4(/) 


(C- 13) 


The next step is to find the moment of z 2 . Equation (C-l) gives 

M M M N 2 M M 


/ ivi jvi m s. z m m M A / 

= (E r ~ + r >i - 6 E r2 a r *i ) = E E r « r c> + E E r >i 

V t = l t = l i= 1 ' t=lj=l irz 1 j -l 


M M 


M M 


- 12 EE r ^- r * 2 ; - i2 EE r » 

i=l j = 1 * = 1 j = l i=l j = l i=l j = \ 


2 2 
r v 

' CJ 1 9] 


(C-14) 


Let z 2 — S\ + 52, where S 1 is the sum of all the terms when i = j and £2 is the sum of all the terms for i ^ j. Then for 
* = 3 


4 4 4 4 

r ci r cj = r ,i r ,j = ™S 

(C-l 5a) 

r ti r tj = ™44 

(C-15b) 

r 6 2 _ r 6 „2 _ 
r ei r tj - r si r cj ~ m 62 

(C-l 5c) 


and 


S 1 = 2M(m 8 -f 19m44 - 127T7 6 2) 


(C-16) 



For i ^ j, let k =| i - j |. Note that r 2 ,r 2 ; - = m\ because the random data make the second moments of r 2 - independent 
of r 2 ,. The double sums can be converted into single sums as follows: 


M-l 


M- 1 


M-l 


M-l 


Si = 2 Y, c * r ci r cj +36 C k r ci r cj r s i r j; + 2 X! c *44 - 24 n c *44 r 


ib=rl 


ib=l 


fc=l 


Jt=l 


M-l 

= 4 c k d k + 36(M - l)Mml - 24(M - l)Mm 4 m| 

k= 1 


(C-17) 


where c * — 2(A/ — &) for & = 1 , 2, 
found from the relation 


M-l 

..M-l. Note that £ c k ~ M 2 - M (and c 0 = A/). The variance of * can now be 
k = 1 


4 = *2-(*) 2 = S 1 + S 2 -4 


Using the previous equations, it can be shown that 


a 2 = 2 M 


rn 8 -f 19 m 4 4 — 12ra62 


+ 4M 


M-l 


— V Cjtdjt — M rn\ — 9 mi + 6m 4 r 

M " 


Jk= 1 


(C-18) 


(C- 19) 


where 

4 = 44 

such that 

A: =| i-j | 

= (r - - cos 4 ( r - 2 cos 4 ^j) 
= (r 2 - re- 8 "* + /*/ 4) 

and 

/jt = cos40 x * cos4</>j 


At high-loop SNR, the mean and variance become 

Hz, ideal — — 4 MR 


(C-20) 


(C-21) 


(C-22a) 


and 


7 z , ideal 


= 64 A/ R 8 - + 


R R 2 


11 JL 

+ & + R 4 , 


(C-22b) 


which result in 
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SN Rz, ideal — 




(C-23) 


So far, no assumption has been made regarding the statistics of the phase error process. Assuming that the phase error 
is indeed Gaussian with second-order density, 

,, , v 1 ( R(0)ti-2R{T)Mi i+Rm]\ 

p{4>i , ( Pji T ) « 27rv//? 2 (0) _ R2 (7) 6XP { 2(R?(0)-&(t)) ) 


1 ( <t>l ~ 2C(T)<t>i<j>j + <j>] 

^W^cW) exp V 2 ^ (1 - c2(r)) 


(C-24) 


where R(t) is the correlation function and can be expressed in terms of the normalized correlation function C(t) as 
R(t) — cr^C(r) where 

C(r k ) = fl - - ') exp(-1.25flir-i) (C-25) 


Bi is the one-sided loop bandwidth and a? the closed-loop error variance. 


r k = Tij = T s | i- j |=| U - tj 


(C-26) 


With this Gaussian assumption, the expected value of cos b<f) will be e ? Hence, we can substitute in all the above 
equations cos4<^ with e -8<7 * and cos 8$ with e~ 32<7 *. Moreover, 

cos 40 j cos 40j = e~ l6a * cosh(16(J^C(r A; )), i ^ j (C-27) 
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